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Assuming the flat FRW universe in Einstein's gravity filled with New Variable Modified Chaplygin 
gas (NVMCG) dark energy and dark matter having negligible pressure. In this research work we 
analyze the viability on the basis of recent observation. Hubble parameter H is expressed in terms 
of the observable parameters Ho, f2^ and the model parameters Ao, Bo, Co, m, n, a and the red 
shift parameter z. Here we find a best fitted parameter range of Ao, Bo keeping < a < 1 and using 
Stern data set (12 points) by minimizing the x 2 test at 66%, 90% and 99% confidence levels. Next 
we do the joint analysis with BAO and CMB observations. Again evaluating the distance modulus 
p(z) vs redshift (2) curve obtained in the model NVMCG with dark matter with the best fitted 
value of the parameters and comparing with that derived from the union2 compilation data. 
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I. INTRODUCTION 

Modern observations like high redshift survey of SNe la CMBR 0-1], WMAP [TMl illus- 

trate the accelerating phase of the universe. Accelerating phase phenomenon of the universe is not 
understood by the standard Big Bang model of cosmology giving a time like singularity in the past. 
Cosmological constant A with the equation of state p = —p is the simplest candidate in Einstein's 
gravity which describe the present accelerating phase of the universe. The aspects of inflation and 
the cosmological constant are not understood well. There are various candidates to play the role of 
the dark energy having negative pressure and describe the accenting view of the present observation. 
The mystifying fluid namely dark energy is understood to dominate the 70% of the Universe which 
violets the strong energy condition and 30% dark matter (cold dark matters plus baryons). Various 
effective candidate of dark energy namely Chaplygin gas with equation of state p = — — , B > 

introduced by fl5j . Furthermore it has been generalized to the form p — — Jlr, < a < 1 [lfjl ] 
and modified to p — Ap — [13, [Hj]. The another candidate of dark energy was introduced by 
Chakraborty et al 19], known as New Variable modified Chaplygin Gas (NVMCG) which follows the 
equation p — A(a) — B \f , < a < 1, a being the scale factor which gives interesting physical significance. 

For flat universe having the energy densities for dust like matter and dark energies we need to know 
the value of critical energy density and H(z) at high accuracy [2(|. The MCG best fits with the 3 year 
WMAP and the SDSS data with the choice of parameters A = 0.085 and a = 1.724 [HJ which are 
improved constraints than the previous ones —0.35 < A < 0.025 (22J. Here we have assumed the new 
variable modified Chaplygin gas model in flat FRW cosmology. The joint data analysis of stern data 
set with BAO and CMB have been analyzed for this model. The distance modulus vs redshift has been 
examined for our model via redshift- magnitude observational data from Supernova type la (Union 2) . 
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A. Basic Equations for Einstein Gravity 

The Friedmann- Robertson- Walker (FRW) metric is considered as 



ds 2 = -c 2 dt 2 + a 2 (t) 



dr 2 



1 — fcr 2 



+ r 2 (de 2 + sin 2 9dcj) 2 ) 



(1) 



where a(t) is the scale factor of the universe and r, 8, <f> are the dimensionless comoving co-ordinates, 
k be the curvature parameter of space-time metric and takes the values k = 0, +1, —1 for flat, closed, 
open universe respectively. 

We consider a spatially flat universe (k = 0) with dark energy and dark matter (non-interacting). Thus 
the Einstein equations becomes (choosing 8ttG = c = 1) 



H 2 



(2) 



H=--(p + p) 



(3) 



where p = pr>E + Pm the total energy density of the universe and p = pde + Pm the total pressure. 
Pde, Pde are the energy density and pressure for dark energy respectively and p m , p m that for dark 
matter. For non-interacting fluid conservation equations become 



and 



p DE + 3H(p DE +Pde) = 



3H(p m + p m ) = 



(4) 



(5) 



We have assumed that the universe is filled with New variable modified Chaplygin Gas (NVMCG) as 
dark energy whose EoS is [19( 



p = A(a)p DE - with < a < 1 

Pde 



(G) 



where A(a), and B(a) are function of the scale factor a. In particular, choose A(a) = Aoa~" 
and B(a) = Boa~ m with Ao, Bo, m, n are positive constants. For n = m = 0, this model reduces 
to modified Chaplygin Gas and for n = 0, the model reduces to the variable modified Chaplygin gas model. 

Expression for the energy density for NVMCG model is obtained from (4) as fl9| 



,'3A a- n 

p = a exp 



Co 



Bo / 3A (I + q) 
Ao V n 



3(l + c) + n 



n /m-3(l + a) 3A (l + a) 
1 I , a 



(J) 



where Co is an integration constant and T(s,t) is the upper incomplete gamma function. 



For dark matter, the EoS is p m — and so p m = p® n a 3 , p°„ = 3£7°„iJo- The Hubble parameter (H) 
in terms of redshift parameter (z) can be expressed as (from eq. (2)) 



3H 2 = (1 + zf 



3H n° m + exp 




_ , Bo (3Ao{l + a) 
Ao \ n 



3(l + c) + r, 
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xr 



m-3(l + a) 3A (l + a) { 1 



l + z 



(8) 



Subsequently, we investigate the bound on the model parameter by observational data fitting. The 
parameters are determined by H(z)-z (Stern), BAO and CMB data analysis [2_4]-|28| using % 2 minimization 
technique from Hubble-redshift data set. 

II. OBSERVATIONAL DATA ANALYSIS H(z)-z (STERN), BAO AND CMB DATA AS A 

CONSTRAINING TOOL 

From the above expression we can write, the Hubble parameter H(z) which can be put in the form as 



where, 



H 2 (A ,B ,a,m,n, C ,z) = H„E 2 (A , B ,a,m,n,C , z) 



E(A ,B ,a,m,n,C ,z) 



(1 + zfl 2 



3H Q n° m + exp 



3^o (t^ 



(9) 



C + 



Bo / 3A,(1 + , 
A V n 



H(l + a) + n-r. 



m-3(l + a) 3A (l + a) ( 1 



l + z 



1/2 



(10) 



Now E(Aq, Bq, a, m, n, Co, z) contains six unknown parameters Aq, Bq, Cq, m, n and a. Now we will 
fixing two parameters and by observational data set the relation between the other two parameters will 
obtain and find the bounds of the parameters. 



z Data 


H(z) 


a 


0.00 


73 


± 8.0 


0.10 


69 


± 12.0 


0.17 


83 


± 8.0 


0.27 


77 


± 14.0 


0.40 


95 


± 17.4 


0.48 


90 


± 60.0 


0.88 


97 


± 40.4 


0.90 


117 


± 23.0 


1.30 


168 


± 17.4 


1.43 


177 


± 18.2 


1.53 


140 


± 14.0 


1.75 


202 


± 40.4 



Table 1: The Hubble parameter H(z) and the standard error a{z) for different values of redshift z. 



A. Analysis of H(z)-z (Stern) data set 

For given a, m,n,Co,z, Aq and Bq can be best fitted by minimizing Xh- z given by 

2 / A T3 n \ {H(A ,Bq, a, m, n, Cq, z) — H i, s (z)) 2 



(11) 
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H(z)-z data (Stern) 
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Figs. 1 and 2 show the variation of An with Bo for a = 0.0001 and a = 0.01 respectively for different confidence 
levels. The 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) contours are plotted for the H(z)-z 

(Stern) analysis. 



where H b s is the observed Hubble parameter at redshift z and <j z is the error associated with that 
particular observation (see table 1) and H represents the theoretical values of Hubble parameter calculated 
for our model. Here we use the observed value of Hubble parameter at different redshifts (twelve data 
points) listed in observed Hubble data by Stern et al [23( to analyze our model. We consider the present 
value of Hubble parameter Ho = 72 ±8Kms^ 1 Mpc^ 1 and a fixed prior distribution. By fixing the model 
parameter a £ [0,1], m, n and Co, we determine the range of other two parameters Aq and Bq by 
minimizing (11). The probability distribution function in terms of the parameters Aq, Bq, Cq, m, n and 
a is given by 



L = j e-^ X 2 H- z P(H )dHQ (12) 

where P(Hq) is the prior distribution function for Hq. Now our best fit analysis with Stern observational 
data support the theoretical range of the parameters. In figures 1 and 2, we plot the graphs for different 
confidence levels 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) contours for a — 0.0001 
and 0.01 respectively and by fixing the other parameters. The best fit values of Aq, Bq and the minimum 
values of x 2 are tabulated in Table 2. 



a 


Aq 


Bq 


v 2 . 

Amm 


0.0001 
0.01 


0.0000869 
0.0000868 


2.988 
3.282 


32.164 
31.925 



Table 2: H(z)-z (Stern): The best fit values of Aq, Bq and the minimum values of \ 2 for m = 7, n = 13, 
Co = 0.1 and for different values of a. 



B. Analysis of H(z)-z with BAO Peak Parameter 

In this section we use the method of joint analysis, the Baryon Acoustic Oscillation (BAO) peak 
parameter proposed by Eiscnstcin et al [29j . Acoustic peaks occurred because cosmological perturbations 
excite sound waves in initial relativistic plasma in the early epoch of the Universe. Sloan Digital Sky 
Survey (SDSS) survey is one of the first redshift survey (46748 luminous red galaxies spectroscopic sample) 
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Figs. 3 and 4 show the variation of Aq with Bo for a = 0.0001 and a = 0.01 respectively for different confidence 
levels. The 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) contours are plotted for the H(z)-z 

(Stern) + BAO joint analysis. 



by which the BAO signal has been directly detected at a scale ~ 100 MPc. The corresponding comoving 
scale of the sound horizon shell is about 150 Mpc in radius. We shall investigate the two parameters Aq 
and Bq for our model using the BAO peak joint analysis for low redshift (with range < z < 0.35) using 
standard x 2 distribution. The BAO peak parameter may be defined by 



A 



/a 



£7(zi)V3 



[■zi dz 

W E(z) 

Zl 



2/3 



(13) 



where 



n° m (i + Zl ) 3 E( Zl y 



(14) 



Here, E{z) is the normalized Hubble parameter (i.e., E(z) = E(Ao, Bq, a, m, n, Co, z)) and z\ — 0.35 
is the typical redshift of the SDSS data sample. This quantity can be used even for more general models 
which do not present a large contribution of dark energy at early times [3Cj| . Now the x 2 function for the 
BAO measurement can be written as in the following form 



Xbao 



{A - 0.469) : 
0.017 2 



(15) 



where the value of the parameter A for the flat model (k = 0) of the FRW universe is obtained by 
A = 0.469 ± 0.017 using SDSS data set [29| from luminous red galaxies survey. Now the total joint data 
analysis (Stern+BAO) for the x 2 function defined by 



XTot 



+ Xbao 



(16) 



Now our best fit analysis with Stern+BAO observational data support the theoretical range of the 
parameters. In figures 3 and 4, we plot the graphs for different confidence levels 66% (solid, blue), 90% 
(dashed, red) and 99% (dashed, black) contours for a — 0.0001 and 0.01 respectively and by fixing the 
other parameters. The best fit values of Aq, Bq and the minimum values of \ 2 are tabulated in Table 3. 
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Figs. 5 and 6 show the variation of An. with Bo for a = 0.0001 and a = 0.01 respectively for different confidence 
levels. The 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) contours are plotted for the H(z)-z 

(Stern)+ BAO+ CMB joint analysis. 



a 


A Q 


Bo 


X 2 ■ 


0.0001 
0.01 


0.0000868 
0.0000867 


2.989 
3.283 


32.165 
31.924 



Table 3: H(z)-z (Stern) + BAO: The best fit values of Aq, Bo and the minimum values of x 2 for m = 7, 
n = 13, Co = 0.1 and for different values of a. 

C. Analysis with H(z)-z, BAO Peak Parameter and CMB Shift Parameter 

Another dynamical parameter that is used in recent cosmological tests is the CMB (Cosmic Microwave 
Background) shift parameter which is the useful quantity to characterize the position of the CMB power 
spectrum first peak. The CMB power spectrum first peak is the shift parameter which is given by [3ll . |32j 



dz' 







H(z')/Ho 



(17) 



where Z2 is the value of z at the surface of last scattering. WMAP data gives 1Z = 1.726 ± 0.018 at 
z = 1091.3. For CMB measurement \ 2 function can be defined as 



2 _ (K- 1.726) 2 

XCMB — 



(0.018) 2 

and the total joint data analysis (Stern+BAO+CMB) for the x 2 function defined by 



(18) 



2 _ 2 i 2 i 2 

XTot — Xh-z t Xbao ~t Xcmb 



(19) 



Now our best fit analysis with Stern + BAO + CMB observational data support the theoretical range 
of the parameters. In figures 5 and 6, we plot the graphs for different confidence levels 66% (solid, blue), 
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Fig.7 

Fig. 7 shows the variation of distance modulus fj,{z) vs redshift z for our model (solid line) and the Union2 

sample (dotted points). 



90% (dashed, red) and 99% (dashed, black) contours for a — 0.0001 and 0.01 respectively and by fixing 
the other parameters. The best fit values of Aq, Bq and the minimum values of \ 2 are tabulated in Table 
4. 



a 


^0 


Bo 


v 2 . 


0.0001 
0.01 


0.0000869 
0.0000867 


2.987 
3.283 


32.165 
31.926 



Table 4: H(z)-z (Stern) + BAO + CMB: The best fit values of A , Bq and the minimum values of \ 2 
for m = 7, n = 13, Cq = 0.1 and for different values of a. 



III. REDSHIFT-MAGNITUDE OBSERVATIONS FROM SUPERNOVAE TYPE la 



Recent observation of high redshift survey Supernovae Type la indicates that the universe undergoing 
an accelerating phase and gives an evidence of existence of dark energy. Since 1995, two teams of High-z 
Supernova Search and the Supernova Cosmology Project have discovered several type la supernovas 
at the high redshifts PHI]. The observations directly measure the distance modulus of a Supernovae 
and its redshift z [33j . Recent observational data, including SNe la which consists of the gold sample 
which has a 157 supernova 34 and the another set is a combined data set of a 192 supernova |33j |. 
Here we consider 557 data points and belongs to the Union2 sample 35] which are considered in the next. 



The luminosity distance d^(z) and distance modulus /i(z) for Supernovas are calculated by 

dz' 



d L {z) = (l + z)H 



and 



/i(z) = 5 lo£ 



10 



o H(z) 

d L (z)/H 
IMPc 



25 



(20) 



(21) 



The best fit of distance modulus as a function /i(z) of redshift z for our theoretical model and the 
Supernova Type la Union2 sample are drawn in figure 7 for our best fit values of Aq, Bq with the other 
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chosen parameters. From the curves, we see that the theoretical NVMCG model is in agreement with 
the union2 sample data. 



IV. DISCUSSION 



We proposed here the FRW universe filled with dark matter (perfect fluid with negligible pressure) 
along with new variable Modified Chaplygin gas (NVMCG) which is one of the candidates of dark 
energy. We present the Hubble parameter H in terms of the observable parameters fi^, Hq with the 
redshift z and the other parameters like Aq, Bq, Co, a, m and n. We have chosen the observed values of 
9? m = 0.28 and H = 72 Kms" 1 Mpc" 1 . From Stern data set (12 points), we have obtained the bounds 
of the arbitrary parameters Aq and Bq (Table 2) by minimizing the x 2 test and by fixing the other 
parameters m = 7, n — 13, a — 0.1, 0.0001 and Co = 0.1. In this way, we may found the bounds of 
any two parameters by fixing the remaining parameters. Next due to joint analysis of BAO and CMB 
observations, we have also obtained the best fit values and the bounds of the parameters {Aq, Bo) (Table 
3 and 4) by fixing some other parameters m = 7, n = 13, Co = 0.1 and a — 0.01, 0.0001. The best-fit 
values and bounds of the parameters are obtained by 66%, 90% and 99% confidence levels are shown in 
figures 1-6 for Stern, Stern+BAO and Stern+BAO+CMB analysis. The distance modulus fi(z) against 
redshift z has been drawn in figure 7 for our theoretical model of the NVMCG for the best fit values 
of the parameters and the observed SNe la Union2 data sample. Here we show that our predicted 
theoretical NVMCG model permitted the observational data sets. The observations do in fact severely 
constrain the nature of allowed composition of matter-energy by constraining the range of the values of 
the parameters for a physically viable NVMCG model. 
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